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ABSTRACT. Brownian motions on a metric graph are defined. Their generators are 
characterized as Laplace operators subject to Wentzell boundary at every vertex. 
Conversely, given a set of Wentzell boundary conditions at the vertices of a metric 
graph, a Brownian motion is constructed pathwise on this graph so that its generator 
satisfies the given boundary conditions. 



1. Introduction and Main Results 



Since the groundbreaking works of Bachelier [2], Einstein [21 22 1, and Smolu- 
chowski [ 86 ] Qthe theory of the Brownian movement had been established as a central, 
recurrent theme in mathematics and physics. In the sequel the Brownian phenome- 
non stimulated the development of many important ideas and theories. A complete 
description of the history is beyond the scope of this introduction, but in keywords we 
want to mention the following: The construction of Wiener space [89] and Wiener's 
approach of statistical mechanics and chaos [90], Ito's theory of stochastic integra- 
tion [39] and stochastic differential equations p0| , Levy's analysis of the fine struc- 
ture of Brownian motion and his theory of the Brownian local time 1 6TJ 62 1 , Feyn- 
man's path integral (27} with its new view towards quantum mechanics, Kac' work on 
path integrals )44p3| . Towards the middle of the last century there were the works by 
Feller [24]-[261 and Ito-McKean 1 41 42 1 on Brownian motions on intervals (see also 
below), Gross' abstract Wiener spaces |34|, Nelson's work on functional integration 
and on the relation between quantum and stochastic dynamics [65-68], giving new 



momentum to Euclidean and to constructive quantum field theory, e.g. 1 33 78 79 



81 83] and nonrelativistic quantum physics, e.g. [82]. Further we want to mention 
the asymptotics of Wiener integrals and large deviation theory (16 75 1, the theory 
of Dirichlet forms (l}[32}[80}, the development of the Malliavin (63} [64} and Hida 
calculi [36 37], and Bismuth's approach to the Atiyah-Singer index theorem [8,9]. 
In addition, there were important developments in other fields, such as engineering, 
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'it seems that Schrodinger |77) was the first to introduce the notion of a first passage time, (in German 
Erstpassagezeit), i.e., a special type of stopping time, in the continuous time context of the Brownian 
motion process. It is striking that this article and the parallel work of Smoluchowski 1 87 1 has practially 
gone unnoticed in the physics literature, while being cited by statisticians, e.g., |28||84| . 
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biology or mathematical finance, which were triggered by the theory of Brownian 
motion. 

The present article is directly linked to the above quoted works by Feller and Ito- 
McKean. So we want to sketch these in little more detail. In his pioneering arti- 
cles 1 24 -26 1, Feller raised the problem of characterizing and constructing all Brow- 
nian motions on a finite or on a semi-infinite interval. In the sequel this problem 
stimulated very important research in the field of stochastic processes, and the prob- 
lem of constructing all such Brownian motions found a complete solution in the work 
of Ito and McKean pTIHSJ via the combination of the theory of the local time of 
Brownian motion [62], and the theory of (strong) Markov processes |10 17 - 19 38 1. 
The central result of these investigations is that the most general Brownian motion on 
the half line M + is determined by a generator which is (one half times) the Laplace 
operator on M + with Wentzell boundary conditions at the origin, i.e., linear combi- 
nations of the function value with the values of the first and second derivative at the 
origin (with coefficients satisfying certain restrictions, see below). Ito and McKean 
showed in (4TJ — partly based on the ideas of Feller |24]-|26j — how to construct the 
paths of such motions: The boundary conditions are implemented by a combination 
of reflection at the origin with a slow down and killing, both on the scale of the local 
time at zero. The ideas contained this article became one of the roots of their highly 
influential book J42J. 

In recent years, there has been a growing interest in metric graphs, that is, piece- 
wise linear spaces with singularities formed by the vertices of the graph. Metric 
graphs arise naturally as models in many domains, such as physics, chemistry, com- 
puter science and engineering to mention just a few — we refer the interested reader 
to (59| for a review of such models and for further references. Therefore it is natural 
to extend Feller's problem to metric graphs. Stochastic processes, in particular Brow- 
nian motions and diffusions, on locally one-dimensional structures, notably on graphs 
and networks, have already been studied in a number of articles of which we want to 
mention @|13l|2^|29H3T]|35j|58j in this context. 

In previous articles [54 - 57 1, two of the current authors studied the self-adjointness 
of Laplace operators on metric graphs and discussed their spectra. This allowed a 
discussion of the associated quantum scattering matrices. Further properties of the 
semigroups generated by Laplace operators on metric graphs, including a Selberg- 
type trace formula and the problem whether these semigroups are positivity preserv- 
ing or contractive, have been studied in [ 49][5T| . This is one of the motivations of our 
present study since semigroups with these properties typically show up in Markov 
processes. Below we will return to this point, see remark 3.5 The wave equation on 
metric graphs and its finite propagation speed has been discussed in [53]. For suitable 
Laplacians free quantum fields on metric graphs satisfying the Klein-Gordon equation 
and Einstein causality were constructed in f76) . 

In [52] the authors have constructed the paths of all possible Brownian motions (in 
the sense defined below) on single vertex graphs using the well-known Walsh process 
|88j, (see also [|6] |72jf74"l|8"5] l) as the starting point. Furthermore, the relation to the 
quantum mechanical scattering is discussed in detail there. The latter article provides 
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an essential input for the construction of all possible Brownian motions on a general 



metric graph in the sense of definition 2.1 (see below) which we carry out here. 



The article is organized in the following way. In section[2]we set up our framework 
and prove our main results: Theorem 2.5 characterizes all possible Brownian motions 
(in the sense of definition |2.1|) on a metric graph Q in terms of Wentzell boundary 



conditions at the vertices. Conversely, theorem 2.8 states that for every choice of a 
set of Wentzell boundary conditions at the vertices as described in theorem |2.5| one 



can construct a Brownian motion on Q implementing these conditions. Theorem 2.5 
is proved in section [3] As a preparation of the proof of theorem 2.8 we consider in 
section [4] the situation where one is given two metric graphs Q\, Q2 with Brownian 
motions X2, X2 in the sense of definition 2. 1 thereon. If one joins some of the external 
edges of Q\ and Q2 to form a new metric graph Q, it is shown how to construct the 
paths of a Brownian motion X on Q by appropriately gluing the paths of X\ and X2 
together. Theorem |2.8| is proved in section [5] via the procedure of section [4] and the 
results in [52 1, where the paths of Brownian motions on star graphs are constructed 
with methods similar to those of Feller (24^(26} and Ito-McKean |41||42]]. The article 
is concluded in section [6] by a discussion of the inclusion of tadpoles. Furthermore, 
there two appendices: one with a technical result on the crossover times which is used 
in section [4j the other about Feller semigroups and resolvents. 

Given these results, it would be interesting to see whether known results for spe- 
cial cases of Brownian motion or diffusions on metric graphs can be extended to all 
Feller processes. For example, an arcsine law has been proved in [4] for the case of 
a Walsh process on a single vertex graph, for the case of a general metric graph we 
refer to (7j[i3J (for a discussion of local time distributions see also |12"|). In a similar 
vein: What about occupation times on edges for the case of general (local) boundary 
conditions of the type ( |2.5| > at the vertices? Can one say something about large devi- 
ations as done for example for Brownian motions without killing and more generally 
for conservative diffusion processes in [30]? What form does the Ito formula take 
in the case of a diffusion process on a metric graph with a generator subject to the 
boundary conditions (|2.5|>? 
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2. Main Results 

In the present article we shall only treat finite metric graphs, and consider a metric 
graph (Q, d) as being defined by a finite collection of finite or semi-infinite closed 
intervals, some of their endpoints — the vertices of the graph — being identified. 
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See figure [T] for an example of a simple, typical metric graph. The metric d is then 
defined in the canonical way as the length of a shortest path between two points along 
the edges (formed by the intervals), and the length along each edge is measured with 
the usual metric on the real line. For a formal definition of metric graphs within 




FIGURE 1. A metric graph Q with 5 vertices, 9 external and 11 in- 
ternal edges. 



the context of graph theory we refer the interested reader, e.g., to [54|55|. Within 
that context our definition above means that we identify — as we may without any 
loss of generality — an abstract metric graph with its geometric graph (see, e.g., 
[43 ]). Moreover, in the sequel it will often be convenient and without any danger of 
confusion to identify an edge of a metric graph with the corresponding interval of the 
real line. Edges isomorphic to M + are called external, while those isomorphic to a 
finite interval — that is, those edges connecting two vertices — are called internal. 
The set of vertices of Q is denoted by V, the set of internal edges by T and the set 
of external edges by £. Moreover we set £ = X U £. The combinatorial structure 
of the graph Q is described by a map 6 from C into V U (V x V) which associates 
with every internal edge i an ordered pair (d~(i), d + (i)) G V x V, d~(i) is called 
the initial vertex of i while d + (i) is its terminal vertex. If i G T is isomorphic to the 
interval [a, b] then d~ (i) corresponds to a, while d + (i) corresponds to b. An external 
edge e is mapped by d to d{e) G V which is the vertex to which e is incident, and 
also in this case we call the vertex the initial vertex of e. 

For the definition of a Brownian motion on the metric graph (Q , d) we take a stand- 



point similar to the one of Knight [48 ] for the semi-line or a finite interval: 



Definition 2.1. A Brownian motion on a metric graph (Q, d) is a diffusion process 
(Xt, t G M+) such that when X starts on an edge e of Q then the process X with 
absorption in the vertex, vertices respectively, to which e is incident is equivalent to 
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a standard one dimensional Brownian motion on the interval e with absorption in the 
endpoint(s) of e. 

Remarks 2.2. By saying that X is a diffusion process we mean that X is a normal, 
strong Markov process (in the sense of [11]), a.s. with paths which are cadlag and 
continuous on [0, Q), where £ is the life time of X. We shall always assume that the 
filtration for X satisfies the "usual conditions". With the help of the well-known first 
passage time formula (e.g., [70] or [42]) for the resolvent of X it is not hard to show 
as in [48 ] that every Brownian motion on a metric graph Q is a Feller process. 

The first crucial problem is then to characterize the behavior of the stochastic pro- 
cess when it reaches one of the vertices of the graph Q, or in other words, the char- 
acterization of the boundary conditions at the vertices of the Laplace operator which 
generates the stochastic process. We want to mention in passing that in an L 2 -setting 
all boundary conditions for Laplace operators on Q which make them self- adjoint have 
been characterized in |54J55J. The first main result of the present paper is Feller's 



theorem for metric graphs. In order to state this theorem we have to introduce some 
notation. 

The Banach space of real valued, continuous functions on Q vanishing at infinity, 
equipped with the sup-norm, is denoted by Cq{Q). We let A denote a universal ceme- 
tery point for all stochastic processes considered, and make the usual convention that 
every / G C {G) is extended to Q U {A} by setting /(A) = 0. 

Consider the generator A of X on Cq(Q) with domain T>{A). Define the space 
Cq(G) to consist of those functions / in Cq(Q) which are twice continuously differ- 
entiable in the open interior Q° = Q \ V of Q, and which are such that their second 
derivative /" extends from Q° to a function in Cq(Q). 

The next lemma, which can be proved with the fundamental theorem of calculus 
and the mean value theorem, states some of the properties of functions in Cq(Q). 
C(v) denotes the set of edges incident with v G V. 

Lemma 2.3. Assume that / belongs to Cq(Q), and consider v 6 V, I G Then 
the inward directional derivatives f^(v{), i = 1, 2, off of first and second order at 
v in direction of the edge I exist, and 

(2-1) /'(«,) = 



lim /'(C)' if vis an initial vertex of I, 

•v,£ei° 



lim /'(C)' if v is a terminal vertex of I, 
(2.2) />j)= lim /"(£) 

hold true. Moreover, f (defined on Q°) vanishes at infinity. 

Remark 2.4. If / G Cg(<?) then by definition of Cg(£7), f"(v k ) = /"(«,) for every 
v G V, and all k, I G C(v), and we shall simply write f"(v). On the other hand, in 
general f(v k ) / f(vi) for k I. 
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Let Vjc denote the subset of V x C given by 

V c = {(v,l), v G V and / G C(v)}. 
We shall also write v\ for (i>, /) G Vc- Consider data of the following form 



a = (a v , v £ V) £ [0,l) V 



(2.3) b=(b Vn v l €Vc)€[0,l] v ^ 

c= (c, v G V) G [0,1] V 

subject to the condition 

(2.4) a v + ^ 

b Vl + c v — 1, for every i> G V. 

Je£(u) 

Define a subspace % a ,&,c of Co (£7) as the space of those functions / in Cq(Q) which 
at every vertex v G V satisfy the Wentzell boundary condition 

(2-5) a v f(v) - Kf\vi) + \ c v f"(v) = 0. 

iec(v) 

Now we can state our first main result: 

Theorem 2.5 (Feller's theorem for metric graphs). Let X be a Brownian motion on 
Q, and let A be its generator on Cq{Q) with domain T>(A). Then there are a, b, c as 



in (fO), dZ4]), so that V{A) = H abc . For f G V(A), Af = 1/2 /". 



Remark 2.6. The boundary conditions in l~L a ^ c are local in the sense that only one 



vertex enters each of the conditions (2.5 1. This is a direct consequence of the path 



properties of X, namely of the condition that the only jumps X may have are those 
from Q (actually from a vertex) to the cemetery point. 

Remark 2.7. Boundary conditions with a v = = c v are often called standard bound- 
ary conditions (see e.g. [56]|59j ) giving rise to what is called a skew Brownian mo- 



tion [42], for a recent survey see e.g. |60|. Killing occurs when a v / 0, and |52| 



provides a detailed discussion of the process for single vertex graphs. When a v = 
the process is conservative and has been studied extensively in [ 29)30) . 



Our second main result is converse of theorem [23] namely 



Theorem 2.8. For any choice of the data as in ( 2.3 1, (|2.4|), there is a Brownian motion 



X on the metric graph Q so that its generator A has T-L a ,b,c os its domain. 



In order to rephrase the statements of theorems 2.5 and 2.8 in a concise way, we 
bring in some additional notation. With a slight abuse of language we shall also call 
any quadruple X = A, P, X) a Brownian motion on Q whenever (fi, A, P) is a 
complete probability space, and X = (X t , t € 1) defined thereon is a Brownian 
motion on Q as in definition |2.1| X{Q) denotes the set of all Brownian motions in 
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this sense, subject to the equivalence relation which is defined by equality of all finite 
dimensional distributions. In consider the (compact, convex) ra-simplex 

n+l 

<7 n = |i £ R n+1 , Xi>0,J2 x i = l}- 

1=1 

Let ctq be the simplex a n , with the point (1, 0, • • • ,0) removed 

^ = a»\{(l,0 > -,0)}. 
It is still convex but not closed. Given a fixed but arbitrary ordering of C(v) any triple 
(a v , (b Vl , I G C(v)), c v ) satisfying ( |2.4| ) can be viewed as an element in 1 with 

n(v) = \C(v)\. With N(g) = Ev( n ( v ) + 2 ) set 

£(£)= X fff^CR^. 

Let i be the mapping defined via theorem[23]by associating to every Brownian motion 
X the data (a, b, c) G Since any two Brownian motions on Q, which have 

the same finite dimensional distributions, define the same semigroup, and therefore 
have the same generator, it follows that i maps these to the same data, that is, i can 



be viewed as mapping from X{Q) to £((7). Theorem 2.8 states that l is surjective. 
To see its injectivity, suppose that [Xi], [X2] are different elements in X{Q), where 
i = 1, 2, denotes the equivalence class of a representative Xj. Assume that 



t([Xi]) = t([X2]). By theorem 2.5 the generator At, i = 1, 2, of X« is uniquely 
determined by the data t([Xj]), and therefore we get A\ = A2. It follows, that Xi and 
X2 define the same semigroup, and therefore all their finite dimensional distributions 
coincide, which is a contradiction. Thus we have proved that 1 is a bijection from 
X(Q) onto 2(g): 

Corollary 2.9. The set X(Q) of all Brownian motions on Q is in one-to-one corre- 
spondence with the setY,(Q). 

3. Proof of Theorem 12.51 



The following notation will be useful throughout this article: If £ is a point in 
Q° = Q \ V, then it is in one-to-one correspondence with its local coordinates (l,x), 
where I G C is the edge to which £ belongs, while x is the point corresponding to £ 
in the interval to which I is isomorphic. Then we simply write £ = (l,x). If / is a 
function on the graph Q we shall also denote /(£) by f(l, x) or fi(x). 

We denote by U = (Ut, t G K+) the semigroup generated by a Brownian motion 
X on Q acting on the Banach space B{Q) of bounded measurable functions on Q, 
equipped with the sup-norm, that is, for / G B(C/), 

U t f(0=Ez(f(X t )), t G K+, £ G 

Clearly, [/ is a positivity preserving contraction semigroup. In the sequel we shall no- 
tationally not distinguish between the semigroup U acting on B{Q) and its restriction 
to the subspace Cq{Q) of B{Q). 
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The proof of the following lemma can be taken over with minor modifications from 
the standard literature, e.g., from (48) Chapter 6.1]. Therefore it is omitted here. 

Lemma 3.1. For every Brownian motion X on the metric graph Q, the generator A of 
its semigroup U acting on Cq(Q) has a domain T>(A) contained in Cq{Q). Moreover, 
for every f € V(A), Af = 1/2 /". 



The preceding lemma implies the second statement of theorem 2.5 The proof of 



the first statement of theorem 2.5 has two rather distinct parts, and therefore we split 



it by proving the following two lemmas: 

Lemma 3.2. Suppose that X is a Brownian motion on a metric graph Q, and that 
T>(A) is the domain of the generator A of its semigroup. Then there are a, b, c as 
in (TO), dZ4}, so that V(A) C Ti a ,b,c- 



Lemma 3.3. Suppose that A is the generator of a Brownian motion X on Q with 
domain T>(A) C T~L a ,b,cfar some a, b, c as in (|2.3|), (|2.4[). Then T>(A) = T-L a ,b,c 



Proof of lemma \3?2\ Our proof follows the one in ]48| Chapter 6.1] quite closely — 
actually, it is sufficient to consider a special case of the proof given there. 

We show that for every vertex v G V there are constants a v G [0, 1), b Vl G [0, 1], 
/ G C(v), c v G [0, 1] satisfying ( |2.4[ >, and such that all / in the domain T>{A) of the 



generator satisfy the boundary condition ( |2.5[ ). To this end, we let / G V {A), fix a 
vertex v G V, and compute Af(v). 

Consider the exit time from v, i.e., the stopping time S v = H{Q°), where for any 
subset M C Q, H{M) = H M denotes the hitting time of M. It is well known 
(e.g., |20][48|[7T| ) that because of the strong Markov property of X, S v is under P v 
exponentially distributed with a rate f3 v G [0, +oo]. Consequently we discuss three 
cases: 

Case j3 v = 0: X is absorbed at v, i.e., v is a trap. Thus Utf(v) = f(v) for all 
t > 0. Consequently, Af(v) = 0, and therefore 1/2 f"(v) = 0. Thus / satisfies the 
boundary condition ( |2.5[ ) at v with a v = 0, c v = 1, and b Vl = for all / G C(v). 



Case < j3 v < +oo: In this case the process stays at v P„-a.s. for a strictly positive, 
finite moment of time, i.e., v is exponentially holding. It is well known (cf., e.g., [48] 



p. 154], [71 p. 104, Prop. 3.13]) that then the process has to leave v by a jump, and by 
our assumption of path continuity on [0, i]), the process has to jump to the cemetery A. 
Therefore we get for t > 0, U t f(v) = exp(- (3t) f (v) , and thus Af(v) + (3f(v) = 0, 
and the boundary condition ( |2.5[ ) holds for the choice 

A 

1 + /T ^ 1 + /3' 



(3-1) a v = - ; - , c,, = - ■ - , b Vl = 0, I e £(v). 



Case f3 v = +oo: In this case the X leaves the vertex v immediately, and it begins a 
Brownian excursion into one of the edges incident with the vertex v. In particular, v 



is not a trap. Therefore we may compute Af(v) in Dynkin's form, e.g., [18 p. 140, 
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ff.], (42] p. 99]. For e > let H ve denote the hitting time of the complement of the 
ball BJv) of radius e around v. Then 



<3 " 2) Am = ] % e„(hJ • 

Now 

E v (f(X(H Vi£ ))^ = Y, fi(e)P v (X(H V;€ ) G l)+f(A)P v (X(H Vi£ ) = A) 
iec(v) 

= £ fi(e)P v [X(H v , e )el), 

l<E£(v) 

where the last equality follows from /(A) = 0. Let us denote 

P v (X(H v>£ )£l) P v {X(H v , t ) = A) 
r ^ e ) = F TTr \ > 1 G r A(e) = 77, x , 

K(e) = l + r A (e) + e ]T r,(e). 

«e£(i>) 

The continuity of the paths of X up to the lifetime ( yields 

£ P v {X(H Vi£ ) G Z) + P v {X(H Vj€ ) = A) = 1, 
ze£(u) 

and therefore equation ( |3.2| ) can be rewritten as 

Iim(A/(t;) + r A (e)/(t;)- ]T r f (e)(/i(e) -/(«)))= 0. 

is£(u) 

Since for all e > 0, K(e)~ 1 < 1, it follows that 

&U(e) Af{v) + K(e) f{v) Jr?K{e) 
which by lemma [XT] we may rewrite as 

\^(a v {e)f{v) + \c v {e)f"{v)- £ b Vl (e) fl(e) ~ f{v) 

iec(v) 

where we have introduced the non-negative quantities 

Observe that for every e > 0, 

a v (e) + c v (e)+ ^(e) = 1. 
ts£(u) 

Therefore every sequence (e n , n G N) with e n > and e n | has a subsequence so 
that a v (e), c v (e) and b Vl (e), I G £(?;), converge along this subsequence to numbers 
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a v , c v , and b Vl respectively in [0, 1], and the relation ( |2.4[ ) holds true. It is not hard to 
check that for every / G Cg {Q) 

m - f(v) 

e 

converges with e \. to f{v{), and therefore we obtain that for every vertex v G V, 
f G "C(vl) satisfies the boundary condition (^5 ) with data a, b, c as in (2.3 1, (2.4). □ 

Before we can prove lemma j3.3| we have to introduce some additional formalism. 

We define the subspace Cq' (Q) of functions / in Cq(Q) which are twice contin- 
uously differentiable on Q° , such that /" (as defined on Q°) vanishes at infinity, and 
furthermore for every v G V and all / G C the limit 

lim /"(£) 



exists. Similarly as in the statement of lemma [23] the last limit is equal to the second 
order derivative f"{vi) of / at v in direction of Z. 

For given data a, b, c as in ( |2.3| ), ( |2.4[ ), it will be convenient to consider T-L a ,b,c 
equivalently as being the subspace of Cq ,2 (Q) so that for its elements / at every 
v G V the boundary conditions ( |2.5| > as well as the boundary condition 

(3.3) f"(vi) = f"(v k ), for all I, k G C(v) 

hold true. Relation ( |3.3| > is just another way to express that /" extends continuously 
from g° to g. 

We consider the sets V, £, and X as being ordered in some arbitrary way. With the 
convention that in C the elements of £ come first this induces also an order relation 
on C. 

Suppose that / G C^ 2 (Q). With the given ordering of £ and X we define the 
following column vectors of length \£\ + 2|X|: 

f(Y) = ((/ e (0), eeS), (/i(0), i 61), (fi{pi) i G X))', 
= ((/e(0), e G £), (/f(0), i G X), (-f'i(pi) i G X 

f00 = ((/e (0), e € £), (/f(0), < G X), i G X 

where the superscript "f ' indicates transposition. 

We want to write the boundary conditions ( |2.5[ ), ( |3.3| ) in a compact way. To this 
end we introduce the following order relation on Vc- For v\, v[, G Vc we set H uj, 
if and only if u -< u' or v = v' and Z ■< I' (where for V and C we use the order 
relations introduced above). For / as above set 

, f 



f(v) = (/(«j), ^ g y £ ) £ , 

/'(y) = (/'(^), vi G ^)*, 
/"(F) = (/>,), u, G Vc)*. 
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Then there exists a permutation matrix P so that 

f(V) = Pf(V), f'(V) = Pf'V), 



f"(V) = Pf"(V). 



In particular, P is an orthogonal matrix which has in every row and in every column 
exactly one entry equal to one while all other entries are zero. 

For every v G V we define the following |jC(u)| x |jC(v)| matrices: 



A(v) 



B(v) 



/ a v 


•• 


• °\ 









•• 


■ 






\0 


•• 


• V 






f-bv h 


~K 




— K, \ 

l \C(v)\ 
















V o 













/ 


(1/2 c v 








• 


• °\ 




1 


-1 





• 


■ 







1 


-1 


• 


• 










1 


-1 • 


• 


1 


V 








• 


• -v 





C(v) 



where we have labeled the elements in C(v) in such a way that in the above defined 
ordering we have l\ ~< I2 ~< ■ ■ ■ -< hc(v)\- Observe that C(v) is invertible if and only 
if c v / 0. Define block matrices A, B, and C by 

A = QA(v), B = QB(v), C = ($C(v). 

vev vev v^v 



Then we can write the boundary conditions ( |2.5[ ), ( |3.3| ) simultaneously for all vertices 
as 

(3.4) Af(V) + Bf'(V) + Cf"{V) = 0. 

Consequently the boundary conditions can equivalently be written in the form 

(3-5) A f{V) + Bf'(V) + Cf"{V) = ^ 

with 

(3.6) A = P~ l AP, B = P^BP, C = P^CP. 

We bring in the following two matrix- valued functions on the complex plane 



(3.7) 



Z±(k) = A±kB + k 2 C, kg 
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Lemma 3.4. There exists R > so that for all k G C with \n\ > R the matrices 
Z±(k) are invertible, and there are constants C, p > so that 



(3.8) 



\Z±(k)~ 1 \\ <C\ 



\k\ > R. 



Proof of lemma 3.4 Since we have 

(3.9) Z±{k) = P~ 1 (A±kB + k 2 C)P 

for an orthogonal matrix P, for the proof of the first statement it suffices to show that 
there exists R > such that 

A ± kB + k 2 C 

are invertible for complex k outside of the open ball of radius R. For this in turn it 
suffices to show that for every vertex v G V the matrices 

A(v) ± kB{v) + k 2 C(v) 



fa v 


± nb Vl + k 2 /2 c v 


±nb Vh 


±nb Vl3 


±nb Vl 




\ 2 


-K 2 













K 2 


— K 2 













K 2 


-K 2 


\ 















±nb, 









are invertible for all k G C with |k| > R. An elementary calculation gives 

*2 



det(A(v) ± kB(v) + k 2 C(v)) 



2\\C(v)\-l 



The choices k = ±1 together with condition \2A) show that the polynomial of sec- 
ond order in k in the first factor on the right hand side does not vanish identically. 
Therefore, it is non-zero in the exterior of an open ball with some radius R v > 0. 
Hence, we obtain the first statement for the choice R = max^gy R v . Moreover, from 
the calculation of the determinants above we also get for every v G V and all k G C 
with \k\ > R an estimate of the form 



(3.10) 



|det(I(u) ± kB(v) + k 2 C(v)) | 1 < const. 



Thus, using the co-factor formula for 

(A(v)±kB{v) + k 2 C{v))~ 1 
we find with ( |3.10| ) the estimate 

\\{A{v) ±kB(v) + K 2 C(v)y 1 \\ < C v \k\ p \ 
for some constants C v ,p v > 0. Consequently we get 

\\(A±kB + ^C)' 1 ]] < C\k\ p , I 



\k\ > R, 



k\ > R, 



for some constants C, p > 0, and by ( |3.9| ) we have proved inequality ( |3.8[ ) 
With these preparations we can enter the 



□ 
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Proof of lemma \33\ Let the data a, b, c be given as in ( |2.3[ ), ( |2.4| ). We have to show 
that the inclusion T>{A) C % a ,bc is not strict. Assume to the contrary that the inclu- 
sion T>{A) C T-L a ,b,c is strict. We will derive a contradiction. 

Let R = (R\, A > 0) be the resolvent of A. Then for every A > 0, R\ is a 
bijection from Cq(Q) onto T>(A), that is, R^ 1 is a bijection from T>(A) onto Co(<5). 
For A > consider the linear mapping H\ : f h-> A/ — 1/2 f" from T-L a ^, c to Co(<?)- 
On X>(.A) this mapping coincides with the bijection i?^" 1 from D(.A) onto Co(£). 
Therefore our assumption entails that H\ cannot be injective. Hence for any A > 
there exists /(A) G H a ,b,c, /(A) / 0, with 

(3-11) ff A /(A) = A/(A)-i/"(A) = 0. 

We will show that /(A) G % a ,f>,c satisfying p. 11) can only hold when /(A) = on 
£/. It will be convenient to change the variable A to k = V2X, and there will be no 
danger of confusion that we shall simply write /(«) for /(A) from now on. Then the 
solution of ( [3.1 1) is necessarily of the form given by 

(3.12) / e («,a?) = r e {K)e~ KX e 6 £, x G M+, 

(3.13) / i (K,x) = r l +( K )e ra + rr( K )e^- a: ) i g X, x G [0, Pi ], 




and we want to show that for some k > 0, the boundary conditions ( |2.5| ) and ( |3.3| ) 
entail that r e (/c) = rf{n) = r~{n) = for all e G £, i G T. For k > 0, define a 
column vector of length \E\ + 2|X| by 

r(«) = ((r e (K), e G £), (r+(«), i G X), ( r r(«), i G X)) 4 , 
and introduce the (|£ | + 2|X|) x + 2|X|) matrices 

X±(k) = 

— appropriately modified in case that £ or X is the empty set — with the |X| x |X| 
diagonal matrices 

e Kp = diag(e Kft , iGl). 
Then the boundary conditions ( |2.5[ ), ( |3.3| ) for /(k) read 

(3.14) Z(«)r(K) = 0, 
with 

(3.15) = (A + k 2 C7)X+(k) + kBX^(k). 

Thus, if we can show that for some k > the matrix Z(k) is invertible, the proof 
of the lemma is finished. Note that the matrix-valued function Z is entire in k, and 
therefore so is its determinant. Thus, if can show that k h-» det Z(k) does not vanish 
identically, then it can only vanish on a discrete subset of the complex plane, and for 
k in the complement of this set Z(k) is invertible. Write 

Xfc(«) = 1±6X(k), 
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with 

so that we can write 
with 





5X(k) =10 



,KfJ 



Z(k) = Z+(k)(1 + 6Z(k)), 



SZ(k) = Z+^k)" 1 Z_(«) 5X{k). 
Observe that in case that X = 0, we obtain 5Z(k) = 0, and in this case the invertibility 



of Z(k) for all k with k > R follows from lemma 3.4 Hence we assume from now 
on that X ^ 0. Lemma |3~4| provides us with the bound 

||Z + (k) _1 ;L(k)|| < const. 

for all k G C, \k\ > R, and for some q > 0. On the other hand, we get 

\\SX(k)\\ <e*">, 

for all k < where po = minj g j pi. Therefore, there exists a constant R' > so that 
for all k < —R' we have ||£Z(re)|| < 1, and therefore for such k, Z(k) is invertible, 
i.e., det Z(k) / 0. Hence there also exists k > so that Z(k) is invertible, and the 
proof is finished. □ 

Remark 3.5. The special (and only) choice of boundary conditions in the form c v = 
and b Vl = (1 — a v )/n(v) > for all v also gives rise to a selfadjoint nonpositive 
Laplace operator on L 2 (Q), see theorem 5.1 in [56]. The associated semigroup is 
positivity preserving and agrees on Cq(Q) n L 2 (Q) with the semigroup associated to 
a corresponding Brownian motion on Q. In turn the theorem just quoted also provides 
examples of selfadjoint Laplace operators, whose semigroups are positivity preserv- 
ing but which are not linked to a Brownian motion process in the above way. 



4. Joining Two Metric Graphs 

For what follows it will be convenient to write Q = (V,X, £, d) for a metric graph 
Q in order to make explicit that V is its set of vertices, X its set of internal and £ its 
set of external edges. 5 is the map from L into V U (V x V) as defined in section[T] 
For simplicity we assume from now on — with the exception of the discussion in 
section|6] — that the metric graphs under consideration do not have any tadpoles, that 
is, internal edges i so that d~{i) = d + {i). 

Throughout this section we suppose that Qk = 04>2fc,£fc, dk), k = 1, 2, are two 
finite metric graphs. In the following subsection we shall construct a new metric graph 
Q = (V, X, £, d) from Q\ and Q2 by connecting some of their external edges. 

It will be convenient to consider the metric graphs Q\, Q2 as subgraphs of the 
metric graph Qq = Q\ l±) Q2 which is their (disjoint) union: Qq = (Vq, Xq, £q, do), with 
Vq = V\ U V%, Xq =X\UX2,£q = £\Ll£2, and where the map do comprises the maps 
d\, 82 in the obvious way. 
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4.1. Construction of the graph Q. Suppose that N is a natural number such that 
N < min(|£i|, l^l)- F° r k = 1, 2, select subsets £1 C of edges with = 
\&2 1 = JV to be joined. Let these sets be labeled as follows 

= {ei, • • • ,ejv}, £ 2 = {'li • • ■ >In}- 




FIGURE 2. Two metric graphs Q%, Q2, to be joined by connecting the 
pairs of external lines (ei, h), (e2, h) and (e3, Z3). 



In addition we assume that we are given strictly positive numbers b±, bjy, 
which will serve as the lengths of the new internal edges, as well as a k £ {—1,1}, 
k = 1, . . . , N, which will determine their orientations. For every k 6 {1, . . . , N} we 
associate with the interval [0, bk] an abstract edge i k (not in Xq) which is isomorphic 
to [0,b k ]. SetX c = {ii, . . . ,i N }, and 

y = Vb, 

I = IoUl c , 

f = £b\(^U^). 

The map d is constructed in two steps: Let & be the restriction of do to 2q U £0 \ (^1 ^ 
£2)- Then 9 is the extension of & to X U 5, which is defined by 

^ fc ) = (! 9i(efc) ' 92( ^' ifafc=1 ' fc=i,...,jv. 

[(52(Zfe),9i(e fc )), if o-fc = — 1, 

Figure[T]shows an example of a metric graph which is constructed from the two metric 
graphs in figure [2]by joining the N = 3 pairs of external edges (ei, (e2,h) and 
(e3, I3). The new internal edges i\, %i and 23 have the lengths 1, y/2 and 1 respectively 
(in some scale). 

Conversely, let a metric graph Q be given. Associate with every vertex v G V of <5 
a single vertex graph with vertex v and n(u) external edges, where n(v) is the 
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number of edges incident with v in Q. Then it is clear that we can reconstruct Q from 
the single vertex graphs Q{v), v G V, by finitely many applications of the joining 
procedure described above. 

For the purposes below it will be convenient to introduce some additional notation. 
We let V c C V denote the subset of vertices of Q which are connected to each other 
by the new internal edges in X c . That is, v G V c is such that there exists at least one 
i G T c with v G d(i). For notational simplicity, here and below we also use d(l) to 
denote the set consisting of d~(l) and d + (l) if I G X, and of d(l) if I G £. In the 
example of the figures [T] and |2j V c = {v2, v$, w\, w^}. 

Consider a vertex v G V c which belongs to Q\, and let ij, G X c , G {1, . . . , iV}, be 
an internal edge connecting v to C/2> i-e., v G d(ik). Then the point 7/ G ^2 w i tn l° ca l 
coordinates (1^, b^) is called a shadow vertex of the vertex shad(w) C Q\ is the 
set of all shadow vertices of v. If v G V C C\Q2, its set of shadow vertices (which are 
points in Q°) is defined analogously. V s = shad(y c ) = U„ e y c shad(f ) is the set of all 
shadow vertices. If £ £ V s , then there exists a unique v G V c so that £ G shad(w). We 
put = w and thereby define a mapping from V s onto V c . Of course, in general 
k is not injective. In figure [3] the shadow vertices of the example above are depicted 
as small circles on the external edges, i.e., V s = {£1, £2, £3, ?7i , ?/2, For example, 
shad(w 2 ) = {»/i,r/2}, shad(u;i) = {£1}, and k(£ 2 ) = ^2, nfa) = V2- 




FIGURE 3. The graphs Q\, Q2 with the shadow vertices (small circles). 



4.2. Construction of a Preliminary Version of the Brownian Motion. From now 
we suppose that we are given a family of probability spaces 

(3°,C°,Q°), £GS , 



and that thereon a Brownian motion with state space Qq in the sense of definition 2.1 



is defined. This Brownian motion is denoted by Z = (Z°(t), t G K+). Actually, 
since Qq = Q\ U £2 and C/i, ^2 are disconnected, this is the same as saying that we 
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are given a Brownian motion on Q% and one on Q2. However, notationally it will be 
more convenient to view this as one stochastic process. We assume, as we may, that 
Z° has exclusively cadlag paths which are continuous up to the lifetime £° of Z°. 
F° = ( J"0, t G R+) denotes the natural filtration of Z°. The hitting time of V s by Z° 
is denoted by r°, i.e., 

r° = inf {t > 0, Z°(t) G V s }. 
Furthermore, we assume that $ = t G M + ) is a family of shift operators for Z° 
acting on S°. 

For any topological space (T, T) denote by Ca (R+ , T) the space of mappings u 
from R + into TU{A} which are right continuous, have left limits in T, are continuous 
up to the lifetime 

C = inf{* > 0, u(t) = A}, 
and which are such that ui(t) = A implies co(s) = A for all s > t. In particular and 
in the present context, to G Ca(M+, Go) is either continuous from JR + into Qq or it 
has a jump from Q\ or Q2 to A, but there can be no jump from Q\ to Q2 or vice versa. 

We shall make use of some special versions of the process Z°, which we introduce 
now. For every v G V c , Z\ = t G M+) denotes a Brownian motion on Qo 

defined on another probability space (Hj, C\, fil) such that under /xj, Z\ is equivalent 
to Z° under We suppose that Z\ exclusively has paths which start in v and which 
belong to Ca(M + , C?o)- (F° r example, one can use a standard path space construction 
to obtain such a version from C°, Q®, Z ).) The hitting time of V s by Z\ is 
denoted by t^, its lifetime by C^. 

The idea to define the preliminary version Y = (Y(t), t G K+) of the Brownian 
motion on is to construct its paths as follows. Let £ G Q be a given starting point. 
<5 (viewed as a set) has the following decomposition (cf. figure [4]): 

Q = Gi*sg 2 , 

with 

9l = 9i\(elU'--\Je° N ), 

g2 = (G2\{l° 1 u---ul° N ))u(i 1 u...ui° N ). 

Thus we may consider £ instead as a point in Q\ tt) Q2 C ^o- 

We pause here for the following remark: Of course, the convention we make that 
all new open inner edges i\ , . . . , i° N are attached to Q2 is somewhat arbitrary. Just 
as well any subset of them could have been attached to Q\ instead. Even though 
different conventions lead to processes with different paths, the main result of this 
section, theorem 4.13 remains unchanged. It follows that all resulting processes are 
equivalent to each other. 

Let Y start as Z° in £ G and consider one trajectory. (In order to avoid any 

confusion, let us point out that even though £ G Ok, k = 1, 2, the process Z° moves 
in Qk-) If this trajectory reaches the cemetery point A before hitting the set V s of 
shadow vertices, it is the complete trajectory of Y and it stays forever at the cemetery. 
If the trajectory hits a shadow vertex 77 G V s before its lifetime expires, this piece of 
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Figure 4. The starting points of Y. 

the trajectory of Y ends at the hitting time r°. Set v = k(t)), and let the trajectory of 
Y continue with an (independent) trajectory of Z\ until its lifetime expires or it hits a 
shadow vertex, and so on. Figure |5]explains the idea. 




Figure 5 . The construction of the process Y . 

The construction described above is formalized in the following way. Define 

H x = X C x = (g)C Q l = (g)/4, 

«ev c v& v c vev c 
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and view each of the stochastic processes Z\,v G V c , as well as the random variables 
Ty, (I, as defined on this product space. Let 

(E n 1 C n ,Q n 1 Z n ,T n ,C), n G N, n > 2, 

be an independent sequence of copies of 

(E\c\Q\z\r\e), 

where Z 1 = (Z*, v G V c ) and similarly for t 1 , Next set 

OO CJO oo 

S=X?, C = (g)C", Q € = Q?®((g)Q ft ), £eS. 

™=0 n =0 n=l 

The procedure sketched above of pasting together pieces of the trajectories of the 
various processes Z™ is controlled by a Markov chain (K n , n G N) which moves at 
random times (S n , n G N) in the state space V c U {A}. We set out to construct this 
chain ((S n ,K n ), n G N). Define 5i = r°. On {Si = +00}, i.e., in the case when 
C° < t°, set = A. Otherwise define 

^ 1 = K (zV)). 

Observe that since all processes considered have right continuous paths, they are all 
measurable stochastic processes, and therefore the evaluation of their time argument 
at a random time yields a well-defined random variable. Set S2 = +00 on {Si = 
+00}, while 

S 2 = Si + r Kl 

on {Si < +00}. On {S2 = +00} put K2 = A, and on its complement 

K 2 = k(Z 1 Ki (t} Ci )). 
These construction steps are iterated in the obvious way: The sequence 

({S n , K n ), n G N) 

is inductively defined by S n = +00 and K n = A on {5 n _i = +00}, while 

S n = S n -i + T"^ n _ 1 , 

on {S n -i < +00}. 

Note that by construction K n = A, n G N, if and only if S n = +00, and in that 
case K n i = A, S n > = +00 for all n' > n. Thus (+00, A) is a cemetery state for the 
chain ((S n ,K n ), n G N). 

For example with a Borel-Cantelli argument it is not hard to see (cf. also [50]) that 
there exists a set 3' G C so that for all £ G Q, Q^{r!) = 0, and for all u G S \ 3' 
the sequence (S n (u), n G N) increases to +00 in such a way that for all n G N, 
S n (u) < 5 n +i(w) holds when S n (cj) < +00. 

Now we are ready to construct Y = (Y(t), G K+). Let £ G Q = ^1 1±) Q2 be a given 
starting point, and suppose that t G is given. On 3' set Y(t) = A. On 3 \ 3' there 
is a unique n G No so that t G [5 n , S^+l). w i m m e convention Sq = 0. If t £ [0, Si), 
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define Y(t) = Z°(t). If t G [S n , S n +i) for n G N, then necessarily S n is finite, so 
that K n G T4, and we define 

F(i) = z£ n (t-s n ). 

In addition, we make the convention F(+oo) = A. The natural filtration generated 
by F will be denoted by T Y = (J" t y , t G K+). 

It follows from the construction of F that A is a cemetery state for Y. Indeed, 
suppose that uj G H \ H', and that the trajectory F( ■ , u) reaches the point A at 
a finite time ( y (uj). This implies that there is an n G No such that ( Y (uj) £ 
[S n (oj), S n+ i(uj)). Then S n (uj) is finite, and therefore K n (uj) G V c so that F(- ,ui) 
is equal to Zj£ ( • — S n (uj),uj) on the interval [S n (uj), S n +i(u)), and this trajectory 
reaches A before hitting a shadow vertex. Hence rjf = +00 which entails that 
S n+ i(ui) = +00. Consequently after S u (uj) there are no finite crossover times for 
this trajectory, and therefore F( ■ , uj) stays at A forever. Furthermore note that the 
left limit Y(( Y (uj)-, uj) at ( Y (uj) belongs to V . 

In terms of the stochastic process Y the random times S n , n G N, have the fol- 
lowing description. Suppose that Y starts in £ G Q. Then S\ is the hitting time of 
V c . But if £ G V c , then actually it is the hitting time of V c \ {£}, because it hits a 
vertex in V c which corresponds to the first hitting of a shadow vertex, i.e., a point in 
Go different from £, by Z°. In particular, S\ > 0. Similarly, S n is the hitting time of 
V c \ {K n -\} by Y after time 5 n _i. In appendix |A| it is shown that for every n G N, 
5* n is a stopping time with respect to T Y . 

It follows from its construction that Y is a normal process, that is, for every £, £ G, 
Qt(Y(t = 0) = £) = 1. Furthermore, all paths of F belong to C A (M + ,£). Let 
SV be the hitting time of the set of vertices V of Q by F. Then Sy < Si, because 
V c C V and therefore we find that F( • A SV) is pathwise equal to Z°( • A Sy), 
where S v denotes the hitting time of V by Z°. Suppose that the starting point £ 
belongs to 1° , I G X U £ , and I is isomorphic to the interval /. Then by definition of 



Z° (cf. definition 2.1 1, the stopped process Z°( ■ A Sy) is equivalent to a standard 
Brownian motion on the interval I with absorption at the endpoint(s) of /. Hence the 
same is true for F: F( • A Sy) is equivalent to a standard Brownian motion on / with 
absorption at the endpoint(s) of /. 

4.3. Markov property of F. For any measurable space (M, M), B(M) denotes 
the space of bounded, measurable functions on M. Every / G B(Q n ), n G N, is 
extended to (Q U {A})™ by ...,&) = 0, ...,£„)€ (0 U {A}) n , whenever 
there is an index k G {1, . . . , n} so that ^ = A. In this subsection we shall prove the 
following 

Proposition 4.1. F has the simple Markov property: For all f G B(Q), s, t G E+, 
£GS, 

(4.1) £ 5 (/(F( S + t))|jf) =£ y(s) (/(y(t))) 

/zo/tfa true Q^—a.s. on {Y(s) 7^ A}. 



BROWNIAN MOTIONS ON METRIC GRAPHS 



21 



The proof of proposition 4.1 is somewhat technical and lengthy. Therefore it will 
be broken up into a sequence of lemmas. 

For every n G N the probability space (E,C,Q^), £ G Q, underlying the con- 
struction of the process Y may be written as the product of the probability spaces 
(3 



\C^-\Qf n 


_1 ) and (E^ n ,C~ n , Q~ n ) with 








n-l 




oo 




= X 




= x^ 




3=0 








n-l 




oo 


c <n-l 




C- n 














n-l 




oo 


Qf" 1 


= Q^((g)Q i ), 


Q- n 


= ®^' 




i=i 







Introduce a family ,6 = (Z3„, n G No) of sub-u-algebras of C by setting 

Obviously, the family B forms a filtration. Furthermore, from the construction of K n 
and S n it is easy to see that the chain ((S n , K n ), n G N) is adapted to £>. 

First we study the chain ((S n , K n ), n G N) in more detail. Recall our convention 
that 5o = 0. We set Bq = {0, E}, and under the law Q v , v G V c , we put Kq = v. 
g G B(R + x V c ) is extended to 1+ x (V c U {A}) by g(+oo, • ) = g( ■ , A) = 
s(+oo, A) = 0. For g G B(R+ x V c ), n G No, define 

(4.2) (U n g)(s,v)=E v (g(s + S n ,K n )), seR+,v£V c . 

Note that [7 = id, and that for every g G B(R + x V c ) and all n G N, f/ n 5 G 
£?(R + x V c ). In particular, the convention mentioned above applies to U n g, too. 

Lemma 4.2. 

(a) For all m, n G N, m < n, £ G Q, s > 0, a«cf every 5 G x 14) f/ie 
following formula holds true Q^—a.s. 

(4.3) E 6 (g(s + SW, if n ) | B m j) = (U n - m g)(s + 5 m , K m ). 

(b) (f/„, n G No) forms a semigroup of linear maps on B(R + x V c ). In particular, 
for all (s,v) G R+ x V c under Q v the chain Us + S n ,K n ), n G No) is a 
homogeneous Markov chain with transition kernel 

P((s, v),A)= Q v ((s + S 1 ,K 1 )€A), AG B(R + x V c ). 



Proof. For m = n formula ( |4.3| > is trivial. Consider the case when n > 2, m = n— 1. 
Let A G v G V^, and put A„ = An {ifn-i = G From the 
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construction of S n and K n 

E ( (g(s + S n ,K n );A v ) 



)<n-2 



j w<n 2 1a, (l >n i g{s + 5 n _x + t%-\ ^rVr 1 ))) dQ^- 1 ) dQf 

I Ujf g(s + u + T ,K{Z (T°)))dQi) dQ\ 

J = <n-2 \J S / U=S n -l 

EAE Kn ^ (g(s + u + S u K{j) | u=5n i ; A t 



= E^(U 1 g)(s + S n - 1 ,K n ^ 1 );A v 

where in the last step we used definition ( |4.2| ). If in the preceding calculation we 
replace the event {K n _\ = v} by {K n _i = A}, we get zero on both sides because 
K„_i = A implies K n = A (see subsection |4.2[ ). Thus summation over v G V c gives 

Et(g(s + S n , K n );A) = E^((U ig )(s + S n -i, K n „ x ); A), 

and equation ( |4.3| > is proved for the case where n > 2 and m = n — 1. As a conse- 
quence we get 

(£/„<?) M = E v (E v (g(s + S n ,if n ) | 

= ^(([/l 5 )( S + S n _l,K„_l)) 
= ° C/lfif)(s,u). 

Now the general semigroup relation U n+m = U n o U m , m, n G No, follows by an 
application of Fubini's theorem. 

Finally we show formula ( |4.3[ ) in the general case: 

Ez(g(s + S n ,K n )\B m ) 

= E^(g(s + S n ,iQ | £ n _i) | £ m ) 

= ^((C/ l5 )( S + 5 ri _ 1) ^ n _ 1 )|H m ) 

= o • • • o C/ig)(s + S m ,K m ) | # m ) 

where the multiple composition in the last expression involves n — m operators U\. 



The semigroup property of (U n , n E No) implies formula (|4.3|). □ 



It will be useful to introduce some additional notation. For r G N, let denote 
the set of all increasingly ordered r-tuples with entries in R+. If u G W + and s£R 
we set u + s = (-ui + s, . . . , u r + s) G K r . u < s means that Ui < s for all i = 1, 
. . . , r or equivalently u r < s. The relations u > s, -u < s, and u > s are defined 
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analogously. In particular, when s < u, then -u — s G W + . For r, q G N, and u G M+, 
w G with u < u>i , define 

(u,u>) = . . . ,u r ,wi, . . . ,w q ) G 

Furthermore, Y(u) stands for (Y(ui), . . . , Y(u r )), and similarly for Z™(u), n G No, 
v G V c . 

In the sequel we shall consider random variables W m (h,g,u) of the following 
form 

(4.4) W m (h,g,u) = h(Y(u)) Xm (u) g(S m+1 , K m+1 ), 

where m G N, h belongs to B{G r ), r G N, g to £(K+ x V c ), and it G M+. Here we 
have set 

Xm{u) = l{S m <«<S m+ i}- 

For s > with s < -u define 

(4.5) W mjS (/i,5,u) = h(Y(u - s)) Xm(w - s) g(s + S m +i,K m+ i), 
so that W m , s =o(h, g, u) = W m {h, g, u). Moreover set 

(4.6) R m (h,g,u)(s,v) = E v (W m)S (h,g,u)), s G K+, s < u, v G V c . 

For the following it will be convenient to let W m ,s{h, g, u) and R m (h, g,u)(s,v), 
v G V c , be defined for all s G R+. To this end we make the convention that = A 
for all t < 0. Then by W miS (/i, 1, u) = Wm(h, l,u — s) the following formula 

(4.7) # m (/t,l,u)(s + M) = Rm(h,l,u- s)(t,v) 

holds for all m G N, h G £(£ r )» r G N, u G W + , s, t G R+, v G V c . 

Suppose that r, g G N, and that h G -B(£ r )» / G B(G q ). Then /i® / denotes the 
function in B{Q r+q ) given by 

(4 8) h ® fi^---^r+q) 

= h(rn,...,r] r )f(r)r + i,...,r)r+q), (m, ■ ■ ■ ,Vr+q) e £ r+9 - 

Lemma 4.3. Suppose that r G N, w G R+, and /i G B(Q r ). 

(a) /fg G N, to G IR+ w/f/i u r < u>, arcd / G -B(^ <? ), then 
(4.9a) i2o(fc®/,l, («,«;)) = Ro(M(f,w,u r )h, l,u) 

holds true, where M(f, w, s)h G B(Q r ) is given by 

(4.9b) (M(f,w,s)h)( V ) = h(r l )E Tlr (W , s (f,h™)), 1^ r ,0< S < W . 

(b) IfgeB(R + x V c ), then 

(4.10a) Ro(h,g,u) = R (N(g,u r )h, 1, u) 

holds, where N(g, s)h G B(Q r ) is given by 
(4.10b) (7V( 5 , s )/ i )( ?? ) = / t ( ?? )^ r ( 5 ( s + 1 S 1 , J fr 1 )), i)ee r , s>0. 
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Proof. Both statements follow from the Markov property of the Brownian motion Z° 
on Go underlying the construction of Y. We only prove statement (b), the proof of (a) 
is similar and therefore omitted. Using the definition of Rq and the construction of 
Y, we compute for s G M+, v G V c , as follows: 

R (h,g,u)(s,v) = E v (h(Y(u - s)) xoO - s)#(s + -S'l, i^i)) 

= ^(/^(u - S)) l {0 <u- s< r°} 9(8 + T°, k(Z°(t°)))) . 

Recall that J 70 denotes the natural filtration of Z°, and d is a family of shift operators 
for Z°. It follows from the definition of the stopping time t° and the path properties 
of Z°, that on {t° > u T — s} the relation r° = u r — s + r° o i9 Ur _ s holds true. 
Moreover, it is easy to check that on this event we have Z°(r°) = Z°(r°) o t? Mr _ s . 
Therefore 

Ro(h,g,u)(s,v) = E v (h(Z°(u - s)) 1{ <m- s < t °} 

x E v (g(u r + r°, k(Z°(t°))) o tf Ur -s I 



= S„ (h(Z°(u - s)) l{ < u - s <rO} 

x^K- s )(3K+r ,^°(r ))))^ 

= ^(y(u - s)) xo(« - s) E Y {u r - s ) (g(u r + 5i, K{)) 

= Ro{N{g,u r )h, l,u)(s,v), 

and the proof is concluded. □ 

Lemma 4.4. For all m, r € N, h € £(g r ), 5 G x V c ), u G M r +) £ G G, the 

formula 

(4.11) E^(W m (h,g,u)\B m ) =Ro(h,g,u)(S m ,K m ) 

holds Q^—a.s. 

Proof. Observe that both side of equation ( |4.1 1[ ) vanish on the set {K m = A}. Let 

A G B m , v G V c , and set A„ = An {K m = v} G B m . Then 

E e (W m (h,g,u);A v ) 

[ A« ( / h i Y ( U )) l {S m <u<S m+1 } 

x g(S m+1 ,K m+1 ) dQ^ m ) dQf 1 - 1 
U„ (/ h(Z™(u - s)) 1{ <«- S <r^} 

Xg{s + T™,K(Z™(r™)))dQ m ) _ dQf™- 1 



<m-l 



<m— 1 



C<m-1 



A„ ( / - S)) 1{ <„- S <r0} 
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<m-l 



S — Sn 



xg{s + r ,K(Z (T )))dQ° v ) 
= E^E v (h(Y(u-s))xo(u-s)g(s + S 1 ,K 1 )) \ s=s j A 

= E<: (E Km (h(Y(u - s)) X o(u - s) g(s + S u Kx)) \ s=Sm , A 

= E^(R (h,g 1 u)(S m ,K m ); A„). 
Summation over v G V c finishes the proof. □ 

Lemma 4.5. For all m, r G N, u G W + , h G B(G r ), g G B(R + x V c ), 

(4.12) R m (h,g,u)(s,v) = (U m Ro(h,g,u))(s,v), s G K+, f G V c , 

Proof. By definition of C/ m 

(U m R (h,g,u))(s,v) = E v (R (h,g,u)(s + 5 m! iC m )). 
With formula ( |4.7| ) and lemma [441 we find 

(y m R (h,g,u))(s,v) = E v (R (h,g(s + • ),u - s)(S m ,K m )) 
= E v (E v (W m (h,g(s+ ■),u-s)\B m )) 
= E v (W m (h,g{s+ -),u-s)) 
= E v (W mtS (h,g,u)) 

= R m (h,g,u)(s,v). □ 

Corollary 4.6. For all m, n, r G N, u G W + , h G £(£ r )> g G 5(R + x F c ), 

(4.13) U n R m (h 
is valid. 

Proof. By lemma [43] and lemma |4~2j statement (b), we obtain 

U n R m (h,g,u) = U n o U m R (h,g,u) = U n+m R (h, g,u) = R n+m (h,g,u). □ 

Lemma 4.7. For all m, n G N, m < n, r G N, u G W + , h G 5(£ r )> g G x l/ c ), 

£ G £/, the following formula holds true: 

(4.14) £€(W n (/i,s,u)|jB m ) =i2 n _ ro (/i, 5 ,«)(5 m ,iir m ). 
Proof. Apply lemma [4~4| to compute as follows 

# e (W n (/t, 5 ,*u)|£ m ) =^(^(^(^,5,^)1^) 

= E^R Q (h,g,u)(S n ,K n )\B m ) 
= {U n -mRo(h,g,u))(S m ,K m ), 



where we used lemma 4.2 formula (|4.3|>, in the last step. An application of lemma 4.5 



concludes the proof. □ 
With these preparations, we are ready for the 
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Proof of Propositioning Assume that / G B(Q), s, t G R + , and that £ G Since 
(S m , m G N) Q^-a.s. strictly increases to +00, and since S m , m G N, is an J 7 ^- 
stopping time (cf., lemma A. 1 in appendix [A]), it suffices to prove that equation ( 4.1 ) 
holds Q^-a.s. for every m G No on {S m < s < S m+ \,Y(s) 7^ A} G Fj . We fix an 
arbitrary m£% Clearly, the family of random variables of the form 

g(Y(w)) l{o<«,<S ro } W m (h,l,u) 

with r, q G N, u G W + , u r = s, w G M^, h G £(£ r ). and 5 G B{Q q ), generates the 
cr-algebra n {Sm < s < S m+ \}. Therefore it is sufficient to show that 



(4.15) 



Et(g(Y(w)) l {0 < w<Sm } W m (h, 1, u) f(Y(s + t))) 

= E^g(Y(w))l mw<Sm} W m (h,l,u)E Y(s) (f(Y(t))))., 



holds for all r, q G N, u G W + with u r = s, w G R q + , h G £(£ r ), 5 G .B^ 9 ) and 
/ G B{Q). (Since the random variables under the expectation signs of both sides of 
equation ( |4.15| ) vanish on the set {Y(s) = A}, we can henceforth safely ignore the 
condition Y(s) / A.) Expand the left hand side of equation ( 4.15| > as follows: 



(4.16) 



E ( (g(Y(w)) l mw<Sm} W m (h, 1, u) f(Y(s + t))) 

00 

= ho<™<S m } W m (h, l,u) W n (f, 1, s + t)). 



Consider the summand with n = m, which is of the form 

E^(g(Y(w)) l {0 < w<Sm } W m {h ® /, 1, (u, s + *))) 

= E^{g{Y(w)) l {0 < w<Sm } Es(W m (h f, 1, (u, s + t))\ B m )) 
= E i {g{Y(w)) l {0 <„,<5 m } Ro(h ® /, 1, (n, s + t))(S m , K m )) , 

where we made use of formula ( |4.1 1 1. Now we apply statement (a) of lemma|43]with 
the choice q = 1 which yields (recall that u r = s) 

Es(g(Y(w)) l { o< w <s m} W m (h <g> /, 1, (u, s + t))) 

= Es(g{Y(w)) l {0 < w<Sm} Ro(M(f, s + t, s)h, 1, u)(5 m , 

= (</(y(t«)) l {0 <»<5 m } ^(W m (M(/, s + 1, s)h, 1, «) I B m )) 

= Et(g(Y(w)) l {0 < w<Sm} W m (M(f, s + t, s)h, 1, «)) , 

where we used formula ( |4.1 1[ ) again. Combining ( |4.4| ) with ( |4.9b| ) in W m (M(/, s + 
t, s)/i, 1, u) we thus have shown 



(4.17) 



Ez(g(Y(w)) l {0 < w <s m } W m (h, 1, u) W m (/, 1, s + i)) 

W m (/», 1, u) E Y(S) ( f(Y(t)) l {0<t <Si})) ■ 
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Next consider a generic summand with n > mon the right hand side of ( |4.16| >. Then 

E^g(Y(w)) l {0 < w<Sm } W m (h, 1, u) W n (f, 1, a + t)) 

= Es(g(Y(w)) l {0 < w<Sm } W m (h, 1, u) Et(W n (f, l,s + t)\ B m+l )) 

= E$(g(Y(w)) l{ <«,<5 m } W m (h, 1, u) i?„_ m _i(/, 1, a + t)(5 m+ i, iT m+1 )) 

= ^(#(^0)) l{0<w<S m }Wm(^-Rn-m-l(/,l,S + t),u)) 



where we used lemma 4.7 in the second step. Conditioning on B m gives 

Es(g(Y(w)) l {0 < w<Sm} W m (h, 1, u) W„(/, 1, s + t)) 

= E c (g(Y(w)) l { o<„,<s m } Et(W m (h, R n - m -i(f, 1, s + t), «) | 

= ^(sC^M) l{0<w<S m } Ro(h, R n -m-l(f, 1, S + t),U)(S m , K m )) 

= E c (g(Y{w)) l {0 < w <s m } R {N(R n -m-i{f, 1, s + a)/i, 1, u)(5 m , tf m )) . 



Here we used lemmas 4.4 4.3 b, and in the last step also u r = s. Applying lemma 4.4 
we get 

Et {g{Y(w)) l {0 < w<Srn} W m (h, 1, u) W n {f, 1, a + tj) 
= Ez(g(Y(w)) l {0 <«,<s m } 

n—m—l 

(/, l,a + t),a)h, l,u) \ B Tl 
= Ez(g(Y(w)) l {0 < w<Sm } W m {N(Rn- m -i(f, 1, a + t),a)h, 1, u)) 
= l{o<u,<5 m } h(Y(u)) Xm(u) 

n—m—l 

For rj £ Q relation ( |4.7| ) yields 

E v {Rn-m-l(f, l,8 + t)(8 + S 1 ,K 1 ))= E v (l^-m-1 (/, 1, t)(Sl, ^l)) 

= £,(£7„(W n _ ro (/,M)|fli)) 
= £„(W n _ m (/,M)), 



with another application of formula ( |4.14[ ). With the choice t] = Y(s) this relation 
therefore gives 

Ez(g(Y(w)) l { o< w <s m } W m (h, 1, u) W n (f, 1, a + t)) 
(4.18) = !{0<«,<5 m } ^(^(w)) Xm(«) 

X #y(s) (/(^(t)) l{5„_ m <t<5„_ m+ i}) ) ■ 
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Formulae (14. 17b and (14. 1 8b entail 



«x5 m | W m (h, 1, u) 1, s + t) 

n=m 

= P f (<?(!») l {0 < w <5 m} W m (^ 1, u) E Y{S) (f(Y(t)) l {0 <t <Sl })) 

oo 

+ S ^(5(^(^)) 1 {0< i ,<5 m }^(^(^))x m («) 
n=m+l 

x £V( S )(/(X(*)) l{s„_ m <t<s„_ m+1 })) 
= P f (g{Y(w)) l {0 < w<Sm } W m (h, l,u) E Y(S) (/(y(t)))) , 



which proves equation ( |4. 1 5 I. 



□ 



4.4. A Brownian motion on Q and its generator. The stochastic process Y and its 
underlying probability family (S, C, Q), Q = (Q%, £ G £/), are not very convenient 
to work with. Therefore we introduce another version in this subsection. As the 
underlying sample space $7 we choose the path space Ca(M+, Q) of Y endowed with 
the cr-algebra A generated by the cylinder sets of C/±(M. + ,G)- Obviously, Y is a 
measurable mapping from (H,C) into (U,A). For £ G Q let denote the image 
measure of under Y. Set P = (P^, £ G £/). Moreover, let the canonical coordinate 
process on (Q,A) be denoted by X = (X t , t G Clearly, X is a version of Y. 
We set X +oc = A, and denote the natural filtration of X by F = {Ft, t G As 
usual Too stands for a {Ft, t G Whenever it is notationally more convenient we 
shall also write X{t) for X t , t G M+. 

Let ily denote the hitting time of the set V of vertices of Q by X: 

H v = m£{t > 0, X t G V}. 

Suppose that X starts in £ G 1°, I G P U £, and that Z is isomorphic to the interval 



I. Then it follows directly from the discussion at the end of subsection 4.2 that the 
stopped process X{ ■ A Hy) is equivalent to a standard Brownian motion on / with 
absorption in the endpoint(s) of /. The necessary path properties of X being obvious, 
we therefore find that X satisfies all defining properties of a Brownian motion on Q 



(cf. definition 2.1 1, except that we still have to prove its strong Markov property. This 
will be done next. 

Let 6 = {6t, t G R_|_) denote the natural family of shift operators on 0: Ot(oj) = 
uj{t + • ) for to G Q. Thus in particular 9 is a family of shift operators for X. Since 
the simple Markov property is a property of the finite dimensional distributions of a 
stochastic process, and the finite dimensional distributions of X and Y coincide, it 
immediately follows from proposition |4. 1 |that X is a Markov process. Then standard 



monotone class arguments (e.g., |47]|7Tj) give the Markov property in the familiar 
general form: 
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Proposition 4.8. Assume that £ G <5, i G M+, a«<i ?/za? W an J 7 ^— measurable, 
positive or integrable random variable on (Q, A, P). Then 

(4.19) E^Wo9 t \T t ) =E Xt (W) : 
holds true P^—a.s. on {JQ 7^ A}. 

A routine argument based on the path properties of X (similar to, but much easier 
than the one used in the proof of lemma A. 1 in appendix [A]) shows that Hy is an 
^-stopping time. We have the following 

Lemma 4.9. X has the strong Markov property with respect to the hitting time Hy. 
That is, for all £ G Q, t G R+, f G B(G), 

(4.20) E^f{X t+Hv ) I T Hv ) = E XHv (f(X t )) 
holds true P^—a.s. 

Proof. To begin with, observe that since Q = Ca{Q) and X is the canonical coordi- 
nate process, there is a natural family a = (at, t G K+) of stopping operators for X, 
namely at(uj) = u( • At), t G M+. Therefore we get that Ft = o-(X sA t, s G M+) for 
any stopping time T relative to T. Indeed, one can show this along the same lines used 
to prove Galmarino's theorem (e.g., [5 , p. 458], (42j p. 86], gSJ p. 43 ff], (7TJ p. 45]). 
Therefore it is sufficient to prove that for all n G N, s\, . . . , s n G M+, t G M+, £ G £/, 
and all 5 G B(Q n ), f G B(<7), the following formula 

A . . . , X( Sl A fly)) + Hy))) 

(4 ' 21 ' = Et(g{ x (*i A #V), • • • , a fly)) s x(ffv) (/(x(t)))) 

holds. Recall that Sy denotes the hitting time of V by Y. Since is the image of 
under Y, and since Sy = Hy o 1", equation ( 4.21[ ) is equivalent to 

(4.22) £? e (G/(y(t + 5 v ))) =^(G^ n5y) (/(K(t)))), 

where we have set 

G = g{Y( Sl ASv),...,Y(s n ASy)). 

Recall that Si denotes the hitting time of V c C V by F, so that Sy < Si. Y 
is progressively measurable relative to T Y which entails that G is measurable with 
respect to "Pg C J 7 ^ C B\ (see also the corresponding argument in the proof of 
lemma AT]). Using the notation of subsection 4.3 we write 

Et(Gf(Y(t + S v ))) 

= E £ (G f(Y(t + Sy); Sy <t + Sy < S^) 
(4-23) 00 

+ £ E^ (G E ( (f(Y(t + u )) Xn (t + u)\ Bi) \ u=Sv ) . 

n=l 



For the last equality — similarly as in the proof of lemma 4.2 — we made use of 



the product structure of the probability space (5, C, Q^), and the fact that Sy < Si, 
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which entails that Sy only depends on the variable lo° G H°. By lemma 4.7 and 



formula (4.7 1 we get for u < Si, n G N, 



E ( (f(Y(t + u)) X n(t + u) | Bi) = Rn-l(f, l,t + u)(5i, K x ) 

= Rn-l(f, l,t)(Si-«, 

Then 

= ^(G J R n _i(/,l,t)(0,^i); 5v = 5i) 

+ ^(Gi?„_ 1 (/,l,t)(S 1 -^,Ki); <5i) 

= £ C (G £ y(Sv ) (/(^(*)) Xn-i(t)) ; Sy = Si) 

+ Ez{GR^ 1 (f,l,t)(S 1 -S v ,K 1 ); Sy < Si), 

because on Sy = Si, Y(Sy) = Y(Si) = Ki. The second term on the right hand 
side of the last equality only involves the random variables A Sv), Si, Sy, and 
K\. They are all defined in terms of the strong Markov process Z° underlying the 



construction of Y (cf. section 4.2 1. Moreover, on the event {Sy < Si} we get from 
the definition of Si as the hitting time of V c that Si = Sy + Si o $g v . Also, on 
{Sy < Si}, Ki = Ki o -d Sv holds true. On the other hand G is measurable with 
respect to J 7 ^, where J 70 is the natural filtration of Z°. Thus the strong Markov 
property of Z° gives 

Es(GR n -i(f,l,t)(Si-S v ,Ki); Sy < Si) 

= E i :(GE Y{Sv) (B n _i(f,l,t)(Si,Ki)); Sy < Si) 
= Es(GE Y{Sv) (f(Y(t))xn(t)); Sy < Si). 

In the last step we used for rj € Q, 

E r) {R n ^i{f,l,t){S u Ki))=E ri (E r) {W n {fA,t)\Bi)) 

= E v [f(Y(t)) Xn (t)), 



with another application of lemma 4.7 and then we made the choice rj = Y(Sy) 



Similarly, for the first term on the right hand side of equation ( |4.23| > we can use the 
strong Markov property of Z° (together with {t + Sy < Si} = {t < Siofl Sv } n 
{Sy < Si}) to show that 



Ez(Gf(Y(t + Sv));Sv<t + Sy < Si) 

= E^GE Y{Sv) (f(Y(t); t < Si)); Sy < Si). 
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Inserting these results into the right hand side of formula ( |4.23| >, we find 

Et{Gf(Y(t + S v ))) 

= Ez(GE Y{Sv) (f(Y(t)); t < Si); S v < S x ) 

oo 

+ J2 E d GE ns v ){f(Y(t))xn(t)); S V < Si) 

n=l 

oo 

+ J2 E d GE Y(S v ){f(Y(t))Xn(%S V = S 1 ) 
n=0 

= Ez{GE Y(Sv) {f(Y(t)))). 
Thus equation ( 4.22| ) holds and lemma 4.9 is proved. □ 
Proposition 4.10. X is a Feller process. 

Proof. It is well-known, that it is sufficient to prove (i) that the resolvent of X pre- 
serves C (6), and (ii) that for all / G C (G), £ G Q, E ( (f(X t )) converges to /(£) 
as t decreases to 0. (A complete proof can be found in appendix [B]) Statement (ii) 
immediately follows by an application of the dominated convergence theorem and the 
fact that X is a normal process with right continuous paths. 

To prove statement (i) consider the resolvent R = (R\, A > 0) of X, and let 
A > 0. Since X is strongly Markovian with respect to the hitting time Hy of the set 
of vertices V (lemma 4.9 1, we get for £ G Q, f G B(Q) the first passage time formula 



(e.g., (7§ or (42)) 

(4.24) (RxfKO = (R?f)(0 + E^e-^ (R x f)(X Hv )), 

where R D is the Dirichlet resolvent. That is, R D is the resolvent of the process 
X with killing at the moment of reaching a vertex of Q. Recall the equivalence of 
the stopped process X(- A Hy) with the Brownian motion with absorption on the 
corresponding interval / stated at the beginning of this subsection. Then we can give 
explicit expressions for all entities appearing in the first passage time formula (|4.24|). 



Using the well-known formulae for Brownian motions on the real line (see, e.g., |20 
p. 73ft], (42j p. 29f]) we find for f, rj G Q°, 

(4.25a) r^{^, V ) = J2 r xA^^hi,nei} + J2 r ^^ 1 U,vee}, 



with 

(4.25b) r ^.(£. r/ ) • 1 ^7, ^ •■-,,-!-:->/, _ ,-\ 2A 



2A fcei. 

where in local coordinates £ = (i, x), rj = x, y G (0, Oj). In the case of an 

external edge e, we get 

(4.25c) r£ e (£, t?) = -L ( e -^»-fl - e -J*(*+v) 

v 2 A ^ 
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with £ = (e,x), r\ = (e,y), x, y G (0 + oo). Remark that both kernels vanish 
whenever £ or rj converge from the interior of any edge to a vertex to which the edge 
is incident. 

Consider the second term on the right hand side of equation ( |4.24| ). Suppose that 
£ G i°, i G T, and that i is isomorphic to [0, Oj]. Assume furthermore that v\, 
V2 are the vertices in V to which i is incident, and that under this isomorphism v\ 
corresponds to 0, while v% corresponds to ctj. Then we get 

Et:{e- XHv (R x f)(X Hv )) = E^e-^i; H Vl < H V2 ) (Rxf)( Vl ) 

+ E i (e- XH ^, H V2 < H V1 ) (R x f)(v 2 ), 

because X has paths which are continuous up to the lifetime of X, and X cannot be 
killed before reaching a vertex. Here H Vk , k = 1, 2, denotes the hitting time of the 
vertex vu- The expectation values in the last line are those of a standard Brownian 
motion and they are well-known, too (see, e.g., (20} p. 73ff], (42] p. 29f]). Thus for 

£ = (i,x), x G [0,ai\, 

Et{e~ XHv (Rxf)(X Hr )) 

(4.26) = sinh(^A( 0i -x)) {Rxf)[vi) + smh(y/2Xx) r^t^s 

sinh(\/2Aai) sinh(v / 2Aaj) 
Similarly, for £ G e° with local coordinates £ = (e, x), x G (0, +oo) we find 

(4.27) Ede~ XHv (R x f)(X Hv )) = e'^ (R x f)(v), 

where v is the vertex to which e is incident. 

With the formulae ( |4.24| >— ( |4~27] > it is straightforward to check that R\ maps C (Q) 
into itself, and the proof is complete. □ 

Since X has right continuous paths, standard results (see, e.g, fTT] Theorem HI.3.1], 
or [9T| Theorem III. 15.3]) provide the 

Corollary 4.11. X is strongly Markovian. 

Thus we have also proved the 



Corollary 4.12. X is a Brownian motion on Q in the sense of definition 2.1 



It remains to calculate the domain of the generator of X, i.e., the boundary condi- 
tions at the vertices. Let v G V, and assume that X starts in v. Then by construction 
of X and Y, X is equivalent to Z° with start in v up to its first hitting of a shadow 
vertex. That is, X is equivalent to Z Q up to the first time X hits a vertex different from 
v. It follows that if v is absorbing for Z°, then it is so for X, and if v is an exponential 
holding point with jump to A, then it is also so for X with the same exponential rate. 
In particular, v is a trap for X if and only if it is a trap for Z°. If v is not a trap, then 
we can use Dynkin's formula p8] p. 140, ff.] to calculate the boundary condition im- 
plemented by X. Clearly, this gives the same boundary conditions as for Z°, because 
Dynkin's formula only involves an arbitrary small neighborhood of the vertex. (See 
also the corresponding arguments in section|3]) Thus we have proved the following 
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Theorem 4.13. X is a Brownian motion on Q whose generator has a domain char- 
acterized by the same boundary conditions as the generator of Z . 

5. Proof of Theorem 12.81 

Suppose that Q = (V,I, £, d) is a metric graph without tadpoles. Let data a, b, c 
as in ( |2.3| > be given which satisfy equation ( |2.4| ). 

With every v G V we associate a single vertex graph Q (v) consisting of the vertex 



v and |£(u)| external edges. In [52| the authors have shown how the construction 



of Brownian motions on a finite or semi-infinite interval by Feller 1 24 - 26 1 and Ito- 
McKean J4T||42| (cf. also [48 ]) can be extended to the case of single vertex graphs. 
For the convenience of the reader, we quickly sketch the method. If b v = and 
c„ = 1 then this trivially is a collection of |£(u)| many standard Brownian motions 
on the real line with absorption at the origin (corresponding to the vertex v), mapped 
onto the external edges of G(v). If b v = and c v < 1 these Brownian motions are 
killed by a jump to A after holding the processes at the origin for an independent 
exponentially distributed time of rate a/c. For b v ^ one uses a Walsh process 
on Q{v) (see [88|, pi), and builds in a time delay as well as killing, both on the 
scale of the local time at the vertex. With appropriately chosen parameters for these 



two mechanisms, theorem 5.7 in |52| states that the so constructed process X v is a 
Brownian motion on Q(v) such that its generator is the 1/2 times the Laplace operator 
acting on / G Cq (Q(v )) with boundary conditions at the vertex v given by ( |2.5[ >. 



Next we build the graph Q by successively connecting appropriately chosen exter- 



nal edges of the single vertex graphs Q(v), v G V, as in section 4.1 Consider the 
stochastic process X which is successively constructed from the Brownian motions 
X v as in subsections 14.2144.41 Theorem 14.131 states that X is a Brownian motion on 
Q which is such that its generator has a domain which is characterized by the same 
boundary conditions at each vertex v G V as for the single vertex graphs G(v). There- 
fore, X is a Brownian motion on Q as in the statement of theorem [278] whose proof is 
therefore complete. 

6. Discussion of Tadpoles 

Suppose that Q is a metric graph which has one tadpole i t connected to a vertex 
v G V. That is, v is simultaneously the initial and final vertex of if. d{it) = (v, v). 
Figure [6] shows a metric graph with a tadpole attached to the vertex v. Let bt be the 
length of i t . Assume furthermore, that we are given data a, b, c as in equations (|2.3|), 



( |2.4[ ). We want to construct a Brownian motion X on Q the implementing boundary 
conditions corresponding to these data. 

Let Q\ be the metric graph obtained from Q by replacing the tadpole by two external 
edges e\, e2, incident with v. Construct a Brownian motion X\ on Q\ corresponding 
to the data a, b, c as above. 

Consider the real line Rasa single vertex graph Q2 with the origin as the vertex vo, 
and edges l\, I2 which are isomorphic to [0, +00), (—00, 0]. Take a Walsh process X2 
on this graph which with probability 1/2 chooses either edge for the next excursion 
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when at the origin. Then X2 is just a "skew" Brownian motion as in [42 p. 115] 
which actually is not skew. That is, it is equivalent to a standard Brownian motion on 
the real line. 

Now join Q\ and Q2 by connecting the pairs {e\,l\), (e2,h) via two new inter- 
nal edges of length b t /2. Denote the resulting metric graph by Q. Figure [7] shows 
this construction. Construct a Brownian motion X on Q from X\ and X2 as in sec- 




h 
h 




Figure 7. Q constructed from Q\ and £2- 

tion [3] By construction, X is equivalent to a standard Brownian motion in every 
neighborhood of vq which is small enough such that it does not include the vertex v. 
Therefore the additional vertex vq of Q does not yield any non-trivial boundary condi- 
tion. Thus, if we identify the open tadpole edge 2° with the subset of Q2 isomorphic to 
( — 6t/l, 6f/2), then we obtain a Brownian motion X on Q implementing the desired 
boundary conditions. 

Obviously, any (finite) number of tadpoles can be handled in the same way. 

Appendix A. On the Crossover Times S n 

We recall from section [3] that in terms of the process Y the crossover times S n , 
neN, can be described as follows. Let Y start in £ € Q. Then S\ is the hitting time 
of V c \ {£} by Y. In particular, S\ > for all paths of Y. For n > 2, S n is the hitting 
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time after 5 n _i of V c \ {K n _i} by Y. Since by construction Y(£ n _i) = K n _\ 
and the paths of Y are continuous on [0, C), we get S n > S n -i for all paths of Y. 
Therefore 

S n = inf {t > S n -i, Y(t) G V c \ {K„-i}} 
= inf {t > S n -i, Y{t) G V c \ {K n ^}} 

holds. 

In this appendix we prove the following 
Lemma A.l. For every n G N, S n is a stopping time relative to T Y . 
Proof. Set S = 0, and for n G N, 

_(V C \{0, ifn = l, 

[ V c \ { A' n _ i } , otherwise . 

For n G N, r G M+ define 

^n,r = {5 n < r < C} 

Bn,r = f| J i 5 "- 1 ^ n ' d ( y (")>^) < 1/m, r < (}. 

meNueQ, 0<ti<r 

We claim that for all n G N, r G R+, 

(A.l) A n>r = B nr . 

To prove this claim, suppose first that co G A n ^ r . Then S n (u) is finite, and therefore 
the set 

{t > 5 n _iH, Y{t,u) g KM} c [o,cH) 

is non-empty. Thus there exists a sequence (m, i G N) in this set which decreases 
to S n (ui). Since Y(-,u) is continuous on [0, C( w ))> it follows that Y(S n (uj), oj) G 
V n {oS). By assumption, S n (uj) < r < C(w), an d therefore the continuity of Y( - ,uj) 
on [0, C(w)) and 5 n _i(o;) < S n (oj) imply that for every m G N there exists u G 
Qfl [S n -i(u)),r] so that d(Y(u,uj),V n (uj)) < 1/m. Hence a; G -B n ,r- 

As for the converse, suppose now that oj G £? n ,r- Then there exists a sequence 
(u m , m G N) in Q n [5 n _i(o;), r] so that d(Y(u m ,u;), V^(w)) converges to zero as 
m tends to infinity. Since (u m , m G N) is bounded we may assume, by selecting a 
subsequence if necessary, that (u m , m G N) converges to some u G [S n _i(o;), r] as 
m — )• +oo. Thus we find that Y(u, u) G V n (u), and therefore Y( ■ , co) hits V n {ui) in 
the interval [S n -\(uS),r\. Consequently, S n (oj) < r, and hence oj G A„ >r , concluding 
the proof of the claim. 

Next we prove by induction that for every n G N, S n is an J^-stopping time. Let 
n = 1. By ( |A.1| ) for every r > 0, 

^V=f| U {d{Y{u),V l )<l/m}n{r<C}. 

meN«£Q,0<«<r 
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Clearly, {r < = {Y(r) G Q} G T Y . Moreover, since V\ is a deterministic set, 
d( • , Vi) is measurable from Q to R + , and therefore {d(Y(u), V\) < 1/m} G J 7 ^ C 
J 7 ^. Hence Ai jr G J 7 ^. Let t > 0, and write 

{Si <i} = {S!<C<i}U A M 

U {5i<r<c}n{C<t})uA M 

:Q,o<r<t 
U Ai, r n {C < *}) U A ht- 

r£Q,0<r<t 

Therefore {Si < t] G J 7 ^, and hence Si is a stopping time relative to T Y . 

Now suppose that n G N, n > 2, and that S^i-i is an J^-stopping time. We 
show that for all r G M+, A„ ir G T Y . First we remark that since Q is a separable 
metric space, the metric d on Q is a measurable mapping from (Q x £/,£>d <8> B^) 
to (K+, For example, this follows from Theorem 1.1.10 in [69], and the 

continuity of d : Q x Q — > M + when Q x Q is equipped with the product topology. 
Consider K n ~\ = Y(S n -\). Since Y has right continuous paths it is progressively 
measurable relative to T Y (e.g., |7l] Proposition 1.4.8]). Thus by Proposition 1.4.9 
in 1 71 1 it follows that K n -\ is J 7 J^ i -measurable. Consequently on < u}, 

K n _\ is J^-measurable. Equation \KA \ reads 

A n , r = f| |J {S n _i < n} n {d(y(«), F c \ < l/m} n {r < c}. 

meN mSQ, 0<«<r 

It follows that A n}r G J 7 ^, as claimed. But then {S n < t] G J 7 ^ for all t G M+ is 
proved with the same argument as at the end of the discussion of the case n = 1. □ 

Appendix B. Feller Semigroups and Resolvents 

In this appendix we give an account of the Feller property of semigroups and re- 
solvents. The material here seems to be quite well-known, and our presentation of 



it owes very much to [70], most notably the inversion formula for the Laplace trans- 



form, equation ( |B.3| > in connection with lemma B.7 On the other hand, we were not 
able to locate a reference where the results are collected and stated in the form in 
which we employ them in the present paper. This applies in particular to the "mixed" 
forms of the statements (iii)-(vi) in theorem [B.3[ which we find especially convenient 
to use in this article. Therefore we also provide proofs for some of the statements. 

Assume that (E, d) is a locally compact separable metric space with Borel a- 
algebra denoted by 13(E). B(E) denotes the space of bounded measurable real valued 
functions on E, Cq(E) the subspace of continuous functions vanishing at infinity. 
B(E) and Cq(E) are equipped with the sup-norm || • ||. 



The following definition is as in 1 71 1: 



Definition B.l. A Feller semigroup is a family U = (Ut, t > 0) of positive linear 
operators on Cq(E) such that 

(i) Uq = id and \\Ut\\ < 1 for every t > 0; 
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(ii) Ut+s = Uto U s for every pair s, t > 0; 

(iii) lim 40 \\U t f - f\\ = for every f G C (E). 

Analogously we define 

Definition B.2. A Feller resolvent is a family R = (R\, A > 0) of positive linear 
operators on Cq(E) such that 

(i) \\R\\\ < X' 1 for every A > 0; 

(ii) R\ — R^ = (p, — \)R\ o R^for every pair A, p > 0; 

(iii) liniA^ \\XR x f - /|| = for every f G C (E). 

In the sequel we shall focus our attention on semigroups U and resolvents R asso- 
ciated with an E'-valued Markov process, and which are a priori defined on B{E). 
(In our notation, we shall not distinguish between U and R as defined on B(E) and 
their restrictions to Cq{E).) 

Let X = (X t , t > 0) be a Markov process with state space E, and let (P x , x G 
E) denote the associated family of probability measures on some measurable space 
(fi, A) so that P x (Xq = x) = 1. E x { ■ ) denotes the expectation with respect to P x . 
We assume throughout that for every / G B{E) the mapping 

(t,x)^E x (f(X t )) 

is measurable from M + x E into R. The semigroup U and resolvent R associated 
with X act on B(E) as follows. For / G B(E), x G E, t > 0, and A > set 

(B.l) U t f(x) = E x (f(X t )), 

roc 

(B.2) Rxf(x)= / e~ xt U t f(x)dt. 

Jo 

Property (i) of Definitions|B.l|and|B.2|is obviously satisfied. The semigroup property, 



(ii) in Definition B.l follows from the Markov property of X, and this in turn implies 



the resolvent equation, (ii) of Definition B.2 Moreover, it follows also from the 



Markov property of X that the semigroup and the resolvent commute. On the other 
hand, in general neither the property that U or R map Cq(E) into itself, nor the strong 



continuity property (iii) in Definitions B.l B.2 hold true on B(E) or on Cq(E) 



If W is a subspace of B(E) the resolvent equation shows that the image of W 
under R\ is independent of the choice of A > 0, and in the sequel we shall denote 
the image by RW. Furthermore, for simplicity we shall write UCq(E) C Cq(E), if 
U t f G C (E) for alU > 0, / G C (E). 

Theorem B.3. The following statements are equivalent: 

(i) U is Feller. 

(ii) R is Feller. 

(iii) UC (E) C C (E), and for all f G C (E), x G E, lim 40 U t f(x) = f(x). 

(iv) UCo(E) C C (E), and for all f G C {E), x G E, limA^, \R\f{x) = 

/04 

(v) RCo(E) C C (E), and for all f G C (E), x G E, lim 40 U t f(x) = f(x). 
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(vi) RC (E) C C (E), and for all f G C (E), x £ E, lim A -+oo XR\f(x) 



e xt 



Remark B.4. The equivalence of statements (i) and (ii) has been shown in 1 14 No. 81, 
p. 291] based on an application of the Hille-Yosida-theorem. 

We prepare a sequence of lemmas. The first one follows directly from the domi- 
nated convergence theorem: 

Lemma B.5. Assume that for f G B(E), Utf — > f as t 4- 0. Then XR\f — > / as 
X — > +oo. 

Lemma B.6. The semigroup U is strongly continuous on RB(E). 

Proof. If strong continuity at t = has been shown, strong continuity at t > follows 
from the semigroup property of U, and the fact that U and R commute. Therefore it 
is enough to show strong continuity at t = 0. 

Let / G B(E), X > 0, t > 0, and consider for x G E the following computation 

U t Rxf(x) - Rxf(x) 

f'OO f'OO 

= / e- Xs E x (f(X t+s ))ds- / e- Xs E x (f(X s ))ds 
Jo Jo 

/ e- Xs E x (f(X s ))ds- / e- Xs E x (f(X s ))ds 
Jt Jo 

fOO ft 

= (e xt - 1) / e- Xs E x (f(X s )) ds - / e~ Xs E x (f(X s )) ds 
Jt Jo 

where we used Fubini's theorem and the Markov property of X. Thus we get the 
following estimation 

(pod rt 
(e A '-l) J e~ Xs ds + e' Xs ds 

= ^(l-e" At ) H/ll, 

which converges to zero as t decreases to zero. □ 

For A > 0, t > 0, / G B(E), x G E set 

00 f_iyH-i 

(B.3) U x f(x) = V 1 j , nXe nXt R nX f(x). 

^— ' n! 

n=l 

Observe that, because of nA||i? n A/|| < ||/||, the last sum converges in B(E). 
For the proof of the next lemma we refer the reader to (70J p. 477 f] : 

Lemma B.7. For all t > 0, / G RB(E), U t x f converges in B(E) to Utf as X tends 
to infinity. 
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Lemma B.8. IfU t C (E) C C {E) for all t > 0, then R X C {E) C C {E),forall 
A > 0. If R\C (E) C C (E), for some A > 0, and R X C (E) is dense in C (E), 
then U t C (E) C C (E)for all t > 0. 

Proof Assume that U t C (E) C C (E) for all i > 0, let / € C (£), x & E, and 
suppose that (x ra , n £ N) is a sequence converging in (E, d) to x. Then a straightfor- 
ward application of the dominated convergence theorem shows that for every A > 0, 
R\f{x n ) converges to R\f{x). Hence R x f £ C (E). 

Now assume that that R\Cq(E) C Cq{E), for some and therefore for all A > 
0, and that R X C (E) is dense in C (E). Consider / £ RC (E), £ > 0, and for 
A > define U^f as in equation ( |B.3[ ). Because i? n A/ £ Cq(E) and the seri es in 



formula ( B.3 1 converges uniformly in x £ E, we get f7 t / £ Cq(E). By lemma B.7 
we find that U^f converges uniformly to Utf as A — > +oo. Hence Utf £ Cq(E). 
Since RCq(E) is dense in C$(E), Ut is a contraction and Cq{E) is closed, we get 
that U t C (E) C CbCE) for every t > 0. □ 

The following lemma is proved as a part of Theorem 17.4 in [46 ] (cf. also the proof 
of Proposition 2.4 in (7TJ). 

Lemma B.9. Assume that RCq(E) C Cq(E), and that for all x £ E, f £ Cq(E), 
limA->oo ^R\f(x) = /(#)• 77ze?i RCo(E) is dense in Cq{E). 

If for all / £ Cq(E), x E E, Utf(x) converges to f{x) as t decreases to zero, 



then similarly as in the proof of lemma B.5 we get that XR x f(x) converges to f(x 



as A — > +oo. Thus we obtain the following 

Corollary B. 10. Assume that RCq(E) C Cq{E), and that for all x £ E, f £ Co(E), 
lim^o Utf(x) = f(x). Then RCq(E) is dense in Cq(E). 

Now we can come to the 



Proof of theorem B.3 We show first the equivalence of statements (i), (ii), (iv), and 
(vi): 



"(i)=^ (h)" Assume that U is Feller. From lemma B.8 it follows that R x Cq(E) C 
Cq(E), A > 0. Let / £ Cq(E). Since U is strongly continuous on Cq{E), lemma B.5 
implies that XR x f converges to / as A tends to +oo. Hence R is Feller. 

"(ii)=^ (vi)" This is trivial. 



(vi)=> (iv)" By lemma B.9| RCq(E) is dense in C (E), and therefore lemma B.8 



entails that UC (E) C C (E). 



(iv)=>- (i)" By lemmas B^and B.9 RC (E) is dense in C Q (E), and therefore by 



lemma B.6 U is strongly continuous on Cq{E). Thus U is Feller. 

Now we prove the equivalence of (i), (iii), and (v): 
"(i)=> (hi)" This is trivial. 



'(ih)=^ (v)" This follows directly from Lemma B.8 
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"(v)=^ (i)" By corollary |BlO} RC (E) is dense in C (E), hence it follows from 
lemma B.8 that UCq(E) C Cq{E). Furthermore, lemma B.6 implies the strong 
continuity of U on RCq(E), and by density therefore on Cq(E). (i) follows. □ 
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